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Abstract 

A quadratic discrete time probabilistic model, for optimal portfolio se- 
lection in (re-)insurance is studied. For positive values of underwriting 
levels, the expected value of the accumulated result is optimized, un- 
der constraints on its variance and on annual ROE's. Existence of a 
unique solution is proved and a Lagrangian formalism is given. An ef- 
fective method for solving the Euler-Lagrange equations is developed. 
The approximate determination of the multipliers is discussed. This 
basic model is an important building block for more complete models. 
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1 Introduction 

Optimal equity allocation and portfolio selection, for a reinsurance company 
with several portfolios (or subsidiaries) leads, in the presence of constraints 
on non-solvency probabilities, market-shares and ROE's, to highly non-linear 
problems. This general situation is studied in a forthcoming reference 0. In 
certain situations it is possible to replace the non-solvency probabilities by 
stronger quadratic constraints (see 0). This gives a new simplified quadratic 
stochastic optimization problem, which solutions (if any) respect the con- 
straints and are approximate solutions of the original problem. The purpose 
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of the present paper is the resolution of the constraint quadratic optimiza- 
tion problem in its most basic setting: Optimization of the expected utility 
(i.e. the final accumulated result) of one single portfolio, under constraints 
on the variance of the utility, on the annual ROE's and on the sign of the 
underwriting levels, which should be positive. This case covers the most basic 
applications and is an important building block in more general situations 
||. The portfolio in this paper is an extension of Markowitz portfolio || 
to a multiperiod stochastic portfolio, as suggested by (c.f. also [jl|). We 
here construct the portfolio such that future results of contracts written at 
different times are distinguishable, which easily allows to consider different 
maturity times. 

We consider a portfolio rj of iV different types of insurance contracts 
concluded at times 0, . . . , T, where T > 1. The amount of the contract of type 
i, where 1 < i < N, being concluded at time t, where < t < T, is denoted 
rjiif). In other words, r/^t) is the number of unit contracts (e.g. the unit is set 
to one FF) of type i. We suppose that the portfolio 77 does not generate any 
financial flows at and after a certain time T + T. Let (AU)(t + 1, 77) be the 
result of the portfolio rj for the period [t, t + l[, let the utility of 77, at time t, 
U(t,r]) = J2 1<s<t (AU)(s,r]), ([/(O,?]) = 0), be the accumulated result for the 
period [0, t[, and let the final utility of the portfolio rj, U (00, rj) = U(T+T, rj) 
be the accumulated result until no more financial flows are generated. (For 
the precise definition of U see formula fl2.1|) .) For 1 < i < N and < t < T, 
i]i(t) and (AU)(t+l, 77) are random variables. In the case under consideration 
the equity K(t) = K(0) + U(t,rj), where K(0) > is the initial equity at 
t = 0. Let Tt be the events which are possible up to time t. The number of 
unit-contracts i]i(t) to conclude at time t, shall be known with certainty at 
time t. This means that r)i(t) is ^-measurable. We impose that the random 
variable i]i(t) has finite variance, for 1 < i < N and < t < T. We now 
introduce the following constraints on the variable 77, (E is the expectation 
operator): 

• C 3 ) E({AU){t + l,rj)) > c{t)E{K{t)), where c(t) G M + is given (con- 
straint on profitability) 

• C 4 ) E((U (00, 77) — E(U (00, r]))) 2 ) < a 2 , where a 2 > is given ( accept- 
able level of the variance of the final utility) 

• Cq) < rji(t), where 1 < i < N (only positive underwriting levels) 

Let C be the set of portfolios 77, such that 77 satisfies constraints (C 3 ), (C4) 
and (Cq) and such that i]i(t) is J^-measurable and has finite second order 
moment, for 1 < i < N and < t < T. We consider the problem of 



Erik Taflin, AXA 
Version 99.07.22 



2 



Equity Allocation and Portfolio 
Selection in Insurance: A simplified . . . 



optimizing the expected accumulated result, i.e. to find all fj G Co, such 
that 



E(U(oo,rj)) = s\xpE(U(oo,7))). (1.1) 

r?GC 



We establish (Corollary |2.6|) under certain mild conditions, (Hi), (H 2 ) and 
(H%) of §2, on the result processes for the unit-contracts, that the optimiza- 
tion problem (|1.1|) has a unique solution, if Co is non empty. Moreover the 
solution, is derived from a Lagrangian formalism (|2.13|) , and is given by for- 



mula ( |2.17| ). Condition (Hi) says that the final utility (sum of all results) 
of a unit contract, written at time k, is independent of events occurred be- 
fore k. In practice, this is generally not true, among other things, because of 
feed-forward phenomena in the prizing. Condition (H 2 ) is equivalent to the 
statement that no non-trivial linear combinations of final utilities, of con- 
tracts written at a certain time, is a certain random-variable. This can also 
be coined, in more financial terms: a underwriting portfolio t](t), constituted 
at time t, can not be risk-free. Condition (H 3 ) says that the final utility of 
unit contracts, written at different times are independent. The conditions 
(Hi), (H 2 ) and (-££3), which excludes many interesting situations, like cyclic 
markets, have only been chosen for simplicity. They can largely be weak- 
ened, without altering the results of this paper. An important point is that 
no particular distributions (statistical laws) are required. 

To determine, practically, the solution r) given by formula ( |2.17| ), we shall 
give, in Appendix [A], an effective method to calculate the inverse of the lin- 
ear operator representing the quadratic form in constraint (C4). We also 
establish that this operator has a finite spectrum. In Appendix [FJ, the de- 
termination of Lagrange multipliers in the solution ( p. 17] ) is studied and a 
simple approximation method is proposed. The proofs of the mathematical 
results, of § ||, Appendix |A| and Appendix [TJ, are given in Appendix 0. This 
paper, which is self contained, is a formalized version of the report ||. 

Acknowledgement: The author would like to thank Jean-Marie Nessi, 
CEO of AXA-Re, and his collaborators for the many interesting discussions, 
which were the starting point of this work. 



2 The model and main results 

We define the portfolio in a probabilistic context, given by a probability space 
(Q, P, T) and a filtration A = {jF 4 } 4gN , of sub a-algebras of the u-algebra T , 
i.e. To = {Q, 0} and T s C T C T for < s < t. The portfolio is composed by 
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N > 1 types of insurance contracts. The utility U{(t, £'), at if G N of the unit 
contract i, 1 < i < N, concluded at time t G {0, . . . , T} is by definition the 
accumulated result in the time interval [t,t'[ for if > t, and Ui(t,t') = 0, for 
< i' < t. We suppose that Ui(t,t') is ^-measurable and that {u{i, t'))t'>o 
is an element in the spaceQ ^(M^), of processes, with finite moments at all 
orders. Since, for given t G {0, . . . , T}, the process t'))t'>o is -4-adapted, 
it follows that (u(t,t'))t>>o G £(R N ,A). The final utility of the unit contract 
i, u°°{t) = Ui(t, s') (= Ui(t, oo)), when the contract does not generate a flow 
after the time s', s' > 0, is T s > measurable. By hypothesis (see §1) there 
exists here s', such that < s' < T + T. We define the utility U (t, rj) of a 
portfolio rj at time t 6 N, where t]i(s) is J^-measurable for 1 < i < N and 
sg{0,... ,T},byQ 

U(t,ri) = Y,v(k)-u(k,t). (2.1) 

k<t 

The result of r\ for the period [t, t + 1[, introduced in §1, is then given by 
(AU)(t + l,rj) = U(t + 1,7]) - U(t,7]). The final utility U(oo,r]) and the 
equity K(t) are given as in §1. Let £^(R N ,A), 1 < q (resp. £ T (R N ,A)) be 
the subspace of elements £ G £ q {R N ,A) (resp. £(M N , A)), with £(t) = for 
£ > T. In optimization problem (|1 . 1| ) it is imposed that the portfolio has 
finite variance, so it is an element of the Hilbert space 7i = £i,(M. N , A), with 
scalar product given by (r],7]') n = Eo<t<T In(( 1 l(t))(^) ■ (v'ityi^dPiu). 

The set Co is well-defined, although the variance of the final utility U (oo, 77) 
is not finite in general, for the (unit-) utility processes (u(t, t'))t'>o G £(M. N , A). 
In fact the quadratic form 

V ^a( V ) = E((U(^, V )f), (2.2) 

in 7i, has a maximal domain T>(a), since for each 77 G 7i, the stochastic process 
(£/(£, ^))t>o is an element of the space £ P (J8L, A), for 1 < p < 2 (which follows 
directly from Schwarz inequality). The solution of optimization problem ( |1 . 1|) 
is largely based on the study of the quadratic form 

77 h- b( V ) = E((U(oo, V ) - E(U(oo, T]))) 2 ), (2.3) 

in Ti,, with (maximal) domain T>(b) = T>(a). 

x Let 1 < q < 00. Then (Xi) <i G if and only if : fi — > l w is J 7 measurable 

and ((^^^(n.R™) = (Jq \Xi{u)\'^ N dP(uj)) 1 / q < 00 for i > 0, where | |kjv is the norm in 
R N . Let £ q (R N , A) the sub-space of A adapted processes in £ q (R N ). We define £(R N ) = 
n^SUR") and £(R N ,A) = n q >!£ q (R N , A). 

2 The scalar product in R N is x ■ y = J2i<i<N x iVi- 
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We shall introduce another optimization problem, having only piece- vice 
linear constraints and which we will prove, to be equivalent to ( |1 . 1|) , in a 
precise way. Since K(0) > 0, it follows from the formula (see §Ij) for the 
equity K(t) and constraint C 3 that E(U(oa,f))) > 0, if the solution 77 ex- 
ists. Moreover, if a > 1 and constraint C 3 (resp. C 6 ) is satisfied, then it 
is also satisfied with 77 replaced by ar/. Therefore, if the solution fj exists, 
then E((U (00, fj) — E(U(oo,f]))) 2 ) = a 2 . With this observation in mind, we 
introduce the optimization problem, which is to find all fj G C\ such that 

b(fj) = inf b( V ), (2.4) 

where C\ is the set of all 77 G 7i such that the following constraints are 
satisfied: 

• C 3 ) E({AU){t + 1,77)) > c(t)E(K(t)), where c(t) G M + is given (con- 
straint on profitability) 

• C4) E(U(oa,r))) > e, where e > is given (acceptable level of the 
expected final utility) 

• Cg) rji(t) > 0, where 1 < i < N (only positive underwriting levels). 

C\ is a closed convex subset of 7i. In fact, the constraint functions in C3 
and C' A are strongly continuous, so it is the intersection of the three closed 
convex cones defined by C' 3 , C' 4 and C' 6 . As we will see, the advantage of this 
formulation is that it is easy to prove an existence result. 

We also introduce the auxiliary optimization problem, which consists of 
finding all 7} G C 2 such that 

a(fj) = inf a( v ), (2.5) 

where C2 is the set of all 77 G Ti. such that the following constraints are 
satisfied: 

• Cg) E((AU)(t + l,r])) > c{t)E{K{t)), where c(t) G R + is given (con- 
straint on profitability) 

• C'l) E(U (00, T))) = e, where e > is given (acceptable level of the 
expected final utility) 

• Cg) f]i(t) > 0, where 1 < % < N (only positive underwriting levels). 
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Ci is a closed convex subset of Tt. 

To solve the optimization problems (|2.4|) and (|2.5| ), we first establish prop- 
erties of the positive quadratic forms b and a. Let b and do be the restriction 
of b and a, respectively, to the domain V(b ) = £>(a ) = £f(R N ,A), where 
£f(R N ,A) = 7i fl £(M N ,A). The corresponding bilinear form is denoted by 
the same symbol b, a, bo and do, respectively, i.e. 

b(£, 77) = E((U(oo, - E(U(oo, 7?) - £(f/(oc, 77)))), (2.6) 

where £,77 G 22(b), etc. 

Lemma 2.1 b (?7), itnt/i domainV{bo) = £f(R N ,A) (resp. a , roi/i domain 
V(ao) = £f{R N ,A)) is a densely defined closeable symmetric quadratic form 
in 7i. 

We make certain (technical) hypotheses on the claim processes: 

• Hi) u°°(k) is independent of for fceN 

• H2) for fceN the N x N (positive) matrix c(k) with elements Cij(k) = 
E({uf(k) - E(u°°(k))(uf(k) - E(u°°(k))) is strictly positive 

• i!/^) uf(k) and Uj°(l) are independent for k 7^ / 

The first hypothesis implies boundedness of the quadratic forms. 

Lemma 2.2 If the hypothesis (Hi) is satisfied, then T>(b) = T>(a) = 7i and 

the quadratic forms b and a are bounded, i.e. there exists C 2 > such that, 

Kv) < C 2 \\v\\h and 0(77) < C 2 \\r]\\ 2 H , (2.7) 

for 77 EH. 

This lemma and the next crucial result show that (the square root of) each 
one of the quadratic forms b and o is equivalent to the norm in 7i. 

Theorem 2.3 If the hypotheses (Hi), (Hi) and (H3) are satisfied, then the 
quadratic forms b and a are bounded from below by strictly positive numbers 
c\ respectively c 2 ,, where c\> c\> 0, i.e. 

b(r}) > C2||t7||^ and a(i]) > Ci||?7||^, (2.8) 

for 77 GW. 
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Let B (resp. A) be the operator in TC, associated with b (resp. a), (by the 
representation theorem), i.e. 

b(£,r)) = (Z,Br)) H (resp. a(£, V ) = (Z,A V ) H ), (2.9) 

for £ e 7i and i] <E H. B and A are strictly positive, bounded, self-adjoint, 
onto 7i and invertible. The inverse operators B~ l : 7i — > 7i and A -1 : 
— ► 7i are bounded self-adjoint operators. There exist c G R, such that 
< cJ < B < A, where I is the identity operator. It follows from formula 
(|2.9| ), that an explicit expression of A is given by 

(A V )(k) = E(t/(oo,r ? ) M °°(A;)|^) (2.10) 

and that an explicit expression of B is given by 

(B V )(k) = E((U(oo, V ) - J E(t/(oo,r ? ))) M 00 (A;)|^), (2.11) 

for r] eH, where < k < T. 

Theorem 2.4 Let the hypotheses (Hi), (H 2 ) and (H%) be satisfied. If C\ 
(resp. C 2 ) is nonempty, then there exists a unique solution fj G C\ (resp. C 2 ) 
of the optimization problem for b (resp. a) 

b(fj) = inf b(rj) (resp. a(fj) = inf a(r])). (2.12) 

»?eCi r?ee 2 

Lemma |2.2j shows that (the square root of) the function in constraint (C4) is 
strongly continuous in Ti and convex. So C is a closed convex subset of 7i. 

Lemma 2.5 Let the hypotheses (Hi), (H 2 ) and (H3) be satisfied. Solutions 
of optimization problems ( |1 . 1\ ) and ( $.4 ) then have the following properties: 
(i) If fj is a solution of equation ([0[) ; then fj is also a solution of equation 
( $■4 ), with E(U(oo,f))) = e, (ii) If f) is a solution of equation ( \2.4 ), then 17 
is also a solution of equation (\1.I\), with a 2 = E((U(oo,r)) — E(U (00, fj))) 2 )- 

This solves the original problem of §|1]: 

Corollary 2.6 Let hypotheses (Hi), (H 2 ) and (H 3 ) be satisfied. If Cq is 
non-empty, then optimization problem ( \1 . i| ) has unique solution fj e C$. 

To construct the solution of the optimization problem ( |2.12| ), in the case 
of b (resp. a), we shall consider a Lagrangian h\ (resp. ), with 
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multipliers. Let Ao, Ai, . . . \f+T-ii an d A* be real numbers and let v G Tt. 
These are the multipliers. A Lagrangian is defined by 

= \<V)~ E ^t(E((AU)(t + l, V ))-c(t)E(K(t))) (213) 

0<t<T+T-l 

- /i(£(£7 (oo, rj)) - e) - (i/, »7) w , 

where 77 G H, with c = b (resp. a). We now have to find the critical points in 
TC, for fixed A, [i and v and determine the multipliers such that the critical 
point in fact is an element of C\ (resp. C2). The multipliers shall satisfy 

H > 0, /i(E(U (00, rj)) -e) = (resp. // G M), (2.14) 



\ t >0, \ t (E((AU)(t + l,fj))-c(t)E(K(t,fj))) = 0, (2.15) 
for < t < f + T - 1 and 

z^(/c) > 0, (vi(k))(u) = a.e. for uj G suppf)i(k), (2-16) 
< k < T and 1 < i < N. Let C\ (resp. C2) be non-empty. 

Theorem 2.7 There are multipliers, satisfying ( \2.14 ), and 
such that the solution 17, of the optimization problem ( \2.12\ ) is given by the 
unique solution 17 of the equation (Dh\ ^ u )(r]) = (resp. {Dh\ „^){rj) = 0). 
Moreover, 

f) = C-\fim+ X th + v), (2.17) 

0<t<f+T-l 

where C = B (resp. C = A) and where in this formula m G 7i is given by 
the linear functional (m, 77)7^ = E(U(oo,r])) and l t G H is given by the linear 
functional (l t ,v)n = E((AU)(t + l,r]))-c(t)E(U(t,r])), forO < t < f+T-1. 

Uniqueness of the multipliers is assured if, for example, m,l , . . . , lf +T _ x 
are linearly independent and c(t) > 0, for < t < T + T — 1. An explicit 
expression of m is given by m(k) = E(u°°(k)\J-'k)- Hypothesis {Hi) then give 

m{k) = E(u°°(k)), (2.18) 

where < k < T. An explicit expression of l t , < t < T + T — 1, is given by 

l t {k) = E{u{k,t + l)-{l + c{t))u{k,t)\F k ), (2.19) 
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if < k < t and by h(k) = 0, if t < k, where in both cases < k < T. We 
note that, if it is supposed in formula (|2.19| ), that u(k,t) and u(k,t + 1) are 
independent of then l t {k) = E(u(k,t + 1) — (1 + c(t))u(k,t)). So when 
this hypothesis, in addition to hypothesis (Hi), is satisfied, then lt(k) is just 



a vector in W N . 



To determine, practically, the solution fj given by formula ( [2.17| ), the 
operator C" 1 is needed. We shall give, in Appendix |A|, an effective method 
to calculate the inverse operator C~ x of C, reminding that C = A (resp. 
C = B) is defined by formula (|2.10|) (resp. ( p.ll|) ) in the case of a (resp. b). 



A Appendix: The operator C 

In this appendix we study the spectral properties of the operator C, defined 
by (|2TT0| ) (resp. fl2TTTD ) when C = A (resp. C = B). We also give an 



algorithm, which determines C~ l . 

The Hilbert space 7i has a canonical decomposition into a direct sum 
T~L = ®o<k<fT^k, where TCk is the quotient space of TL and the subspace of 
elements rj e Ti., having i](k) = 0. To obtain the corresponding decomposition 
of the operators A and B, we introduce the notation 

A(k,l)rj(l) = E(u 00 (k)(u°°(l)- rj(l))\Fk) (A.1) 

and 

B(k, l) V (l) = E{u°°{k){u™{l) ■ rj(l) - E(u°°(l) ■ v(l)))W, (A.2) 

where rj G H, < k < f and < / < f . Formula (|2~T0|) (resp. ( gH|) ) then 
gives that (Arj) (k) = T,Q<i<f A i k ^ l )v{ 1 ) ( res P- ( B v)(fy = Eo<kt B(jk,l) 
r)(l)). The continuity of the operators A and B in 7i shows that v4(fc, Z) and 
S(A;, I) are continuous operators from Hi to Tifc. 

We also introduce the N x N real symmetric strictly positive matrices 
M a (k) and M b (k), for < < f . Here M^(k) = E(uf{k)uf{k)) and 
M b (fc) have the same elements as the matrix c(k) in hypothesis (-£^2)) i-e. 
Afg(jfe) = E((u?°(k) - E{uf{k)){uf{k) - E(u°°(k))), for < fc < A and 
< / < A. According to hypothesis (i^i), it follows that 

A(k,k)r)(k) = M a {k) V {k), (A.3) 

for < k < T. The spectrum a(M a (k)) (resp. a(M b (k))) of the matrix 
M a (k) (resp. M b (/c)) is a finite set of strictly positive real numbers, according 
to hypothesis (-£^2)- Let A G 1R be in the resolvent set of M a (k). The inverse 
(A(k, k) — A) -1 is then simply obtained by inverting the matrix M a (k) — A. 
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Proposition A.l Let A G R be in the resolvent sets of M a (k) and M b (k). 
Then the inverse (B(k,k) — A) -1 of the operator B(k,k) — A exists and is 
given by 

(B(k, k) - A)-V*0 = (M a (k) - \)-\r,{k) - E( V (k))) 

+(M b (k)-X)- 1 E( V (k)), [ ' j 

for < k < f . 

Let C{k,l) = A(k,l) (resp. B(k,l)), for C = A (resp. B). To solve the 
equation Cr\ = £, where C, E His given, we shall study the system of equations 

C(k, l)r){l) = £(h), < k < T . (A.5) 

0<l<f 

The operator C(k,l) is continuous from Hi to Hk- Although, the operators 
C(k,l) are integral operators in general, the system ( |A.5| ) can be solved in 
a finite number of steps only involving linear algebra in ~R N and evaluation 
of expectation values. To study also the spectral properties of the operator 
C in H, we give this algorithm for a more general equation in H, namely 
Crj - Xr] = f , i.e. 

52(C(k,l)-6 u \)ri(l)=£(k), 0<A;<f, (A.6) 

0<l<f 

when A is in a certain subset, to be determined, of the resolventset E — cr(C), 
of the operator C. 

The first part of the algorithm consists of transforming the system ( |A.6| ) 
into an equivalent lower triangular system. For a given real number A outside 
a certain finite set (of poles), we shall define operators C n (k, I) and elements 
£ n G H, where k, I, n G {0, . . . , T}. The definition will be specified by an finite 
iteration, beginning with n = T and ending with n = 0. Let us introduce the 
N x N real symmetric matrices N*(k) and N b (k) and Df(k)) by 

( N n( k ))ij = E{E{u?{k)\F n )E{uf{k)\F n )) (A.7) 

and 

(JVj(*))« = E((E(u?(k)\F n ) - E{u?{k))){E{uY{k)\F n ) - £(tif (*)))), 

(A.8) 

for < i,j < N, < k < T and < n < T. Symmetric matrices 
Df(k), < k < T and real numbers ff , gf, and df are given by 

Df{k) = M a {k) — A, /r = l, g f = 0, (A.9) 
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and 



d T = E(u°°(T)) ■ ((Df(T)y E(u°°(T))), (A.10) 

for A 6 R - cx(M a (k)). For given integer n, where < n < T, given real 
numbers /„, g n and given N x N real symmetric matrices D n (k), where 
< k < n and where D n {n) is invertible, we define real numbers g n -i, 
d n and N x N real symmetric matrices D n _i(k) : 

d n = E(u°°(n)) ■ ((D n (n))-\E;(u°°(n))), (A.11) 

fn-l = /n(l - drJn), £n-l = 9n + d n fni (A.12) 

and 

ZVi(*0 = D n (k) - d n f 2 n N a n (k), < k < n - 1. (A.13) 



Let cZq be given by formula ( |A.11| ), with n = 0. The elements / n , and o? n of 



the sequences fo, . . . , ff, go, ■ ■ ■ ,9f an d c?o ; • • • ,df are real valued rational 
functions of A in K. Similarly, each matrix element of the matrix D n (n), of 
the sequences Dq(Q), . . . , Df(T), is a real valued rational function of A in IR. 
We can now introduce the linear continuous operators C n (k,l) : Hi — > TCk, 
for k, I, n 6 {0, . . . , f }. We define C T (k, I) by 

C T {k,l) = C(k,l), k^l (A.14) 

and 

C T (k, k) = C(k, k) - XI, (A.15) 

(where / is the identity operator on H k ), for k, I G {0, . . . , T} and A G M. 
For given integer n, where < n < T, and given linear continuous operators 
C n (k, I) :Hi^H k , where k, I G {0, . . . , T} we define operators C"" 1 ^, I) : 
7^^ — > 7^^. If A G R is such that, D n (n) is invertible, then we define 

C n -\k, I) = C n {k, I), n < k < f , < I < f , (A.16) 
C n -\k, /) = 0, < k < n, n < I < f , (A.17) 
C n -\k,l) = f n ^C{k,l), k ^ I, < k < n, < I < n, (A.18) 
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A n -\k, k)ri(k) = D n ^(k)r)(k), < k < n, rj(k) G H k (A.19) 

and 

B n -\k,k)r ] (k) = D n . 1 (k)r ] (k) 

- (1 - g n ^)E{u™{k)){E{u™{k)) ■ E( V (k))), (A.20) 
< k < n, rj(k) e H k . 

If C = A (resp. B), then C n ^(k,k) = A n ~\k,k) (resp. C n -\k,k) = 
B n ' 1 (k,k)), for < k < n. 

According to ( |A.19| ), the operator A n (k, k), in H k , has a bounded inverse, 
if A is such that, D n (k) is invertible: 

(A n (k, k)y l V (k) = (D n (k))- l r,(k), < k < n, V (k) G H k . (A.21) 

As in the case of (B(k, k) — A)" 1 , (see formula ( |A.4| )), an explicit expression 
for the inverse of B n (k, k) can be given. 

Proposition A. 2 Let A be in the resolventset of the two operators M c (k) — 
E n +i<Kfi/r^r(^ for c = a and c = b. Then 

(B n (k, k))" 1 ^) 

= (M a (k)- d r f r N?(k)-\y\ V {k)-E( V (k))) 

n+l<r<f (A. 22) 

+ {M\k)- drfiN^k) - X^EWk)), 

n+l<r<T 

where < k < n < T. 
For reference we note that 

D n (k) = M a (k) — d r f?N r a (k)-\, 

n+l<r<T 

according to ( |C.14|) , of Appendix 

Next we define £ n_1 G H. For given £ G H, let £ T = £. For given integer 
n, < n < T, for given linear continuous operators C n {k,l) : Hi — > Hk, 
k,l G {0, . . . ,T}, where C n (n,n) has a bounded inverse, and for given £ n , 
we define 

C-\k)=C(k), n<k<f. (A.23) 
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and 

C~\k) =C(k) -C n {k,n){C n {n,n)y l C{n), < k < n. (A.24) 

To give spectral properties of the operator C, we introduce notations for 
unions of spectra of certain matrices. We recall that, if X is a linear operator 
in a Hilbert space, then cr(X) denotes here its spectrum. For commodity we 
first introduce cr|, = = 0, the empty set. For < n < T, let 

< = < +1 Ua(Af»- Yl d rfrK(n)) (A.25) 



n+l<r<T 



and let 



a b n = a a a Ua(M\n)~ £ d r f r N b r {n)). (A.26) 

n+l<r<T 

We note that cr£ C cr£_ 1; for c = a and for c = b. 

Proposition A. 3 If X ^ °~n+ii then C T , ... ,C n and £ T , ... , £ n , given by 
( \A.Q )- ftA.24) are sequences of linear continuous operators in 7i and elements 
ofTi, respectively and the system ( \A.dj ) has the same solutions rj G as the 
system 

£ C n (k, l) V (l) = C(k), < k < f . (A.27) 

0<1<T 

For n = 0, the system ( |A.27| ) is lower triangular and reads 

J2 C°(k, l)rj(l) = £°(A;), < k < f , (A.28) 

0<l<k 

where we remind that C° and £° are dependent of A ^ o~{. 

The second part of the algorithm consists of solving the system ( |A.28| ) 
by successive substitutions. 

Proposition A. 4 If \ £ (resp. A ^ a^), when C = A (resp. C = B), 
then r](0), . . . ,r](T), are successively given by 

r ] (k) = (C k (k,k))-\e(k)- Yl C (k,l) V (l)), 0<k<T, (A.29) 

0<Kfc-l 

where we note that the sum is absent for k = 0. 

We remind that the operator (C k (k, fc)) -1 , < k < T, is explicitly given by 
(pOTl ) (resp. (pQ2D ), when C = A (resp. C = B). 

To sum up, if A ^ (Tq (resp. A ^ Uq), then 77, given by formula ( |A.29|) , 
is the unique solution, in H, of equation ( |A.6| ), in the case of C = A (resp. 
C = B). Consequently c(A) C (Tq an d o~(B) C Uq. In particular the spectrum 
o~(C) is a finite set of strictly positive real numbers. 
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B The Lagrange multipliers 



We have already established that there exists multipliers fi, A , . . . , Xf+T-i G 
R and v G TC, such that the unique solution fj of the optimization problem 
( 2.12]) is given by formula ( 2.17|) . In this appendix, we derive equations of 



which the multipliers are solutions. These equations are taken as a starting 
point for the derivation of algorithms permitting to calculate approxima- 
tions of the multipliers. Starting with the zeroth approximation as being 
the value of the multipliers corresponding to the solution of the determinis- 
tic approximation of the optimization problem Q2.12 ), we derive an explicit 



first approximation of the multipliers in the stochastic case. More precise 
approximations is a subject of future studies. 

A basic building block in our equations for the multipliers, is the solution 
of a certain finite dimensional optimization problem. (The dimension will 
mostly be N or T + T, i.e. the number of types of contracts or number of 
constraints in (C 3 ) and (C 4 ) of § Q respectively). For z G M n , let z > be 
defined by Z{ > for 1 < i < n. Let m be a strictly positive symmetric 
operator in W 1 . For given x G M. n , let F^(x) > be the unique solution of 
the optimization problem 

1 1 

-F+(x) ■ {mF+(x)) - x ■ F+(x)) = inf (-y ■ {my) -x-y). (B.l) 

Z y>0 Z 

The solution F+(x) of ( |B.1| ) is also the unique critical point of the Lagrangian 
L x '■ 

L x (y) = ^y-(my)-x-y-F~(x)-y, (B.2) 

where 1 < i < n, are the multipliers corresponding to the constraint 

y > 0. We have 



where 



m-\x + F-(x)), (B.3) 



F+(a;)>0, F~(x)>0, F+(x) ■ F~{x) = 0. (B.4) 



We note that there exists, for given m, a function x t— > P x of M n to linear 
projections on K n , (i.e. P x = Pj), taking a finite number (2 n ) of values, such 
that F*(x) = m~ x P x x. 

We shall here consider the solution fj (given by (|2.17|) ) only in the case 
of C = B. The case C = A, is similar but simpler. Let Xf+T = ^ an d let 
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If+T — m - We define Q e (k) by 

© e (*) = *W E M(A;) + (M«(A;)-M t '(A;)) J E;(r)(A;)))-^i?(fc,/)77(/)), 

0<t<T+T l+k 

(B.5) 

for < k <T and e = ±. It follows that 

fj(k) = Q + (k) and u(k) = &~{k), (B.6) 

for < fc < f . 

For simplicity we suppose in the sequel of this paragraph that the subset 
{E(Iq),... ,E(lf +T )} of Ti. is linearly independent. This ensures that the 
multipliers are unique in the deterministic case as well as in the stochastic 
case. Expressions ( P-5|) and (p.6| ) are satisfied by the multipliers and the 



solution f) given by formula ( |2.17| ). To obtain another such equation let L 
be the strictly positive symmetric operator in M T+T defined by (L -1 )^ = 
(l t ,B-H s ) and let r(v) E R T+T be given by (r(v)) t = (h, B~ x p). Then the 
multiplier A satisfies: 

X = F-(L(r(u)-e)). (B.7) 

Equations ( |B.5| ), ( |B.6| ) and ( |B.7| ) form a closed system for the solution f) 
and the multipliers A and v, (reminding that Xf+r = /-O- Alternatively, 
expression ( ]2.17|) of rj, formula ( |B.5| ) with e = — , the expression of v in ( |B.6[ ) 
and expression (p. 7] ) of A, also form a closed system of equations, for rj, X 
and v. 

We shall next give a crude approximation of the multipliers A and v. To 
construct the zeroth approximation, we consider the deterministic approxi- 
mation of the optimization problem ( |2.12| ). Let Cd be the subset of elements 



7] 6 C±, such that rj(k) is JF measurable for < k < T. The deterministic 
approximation, of the optimization problem, is then to find all t)o G Co such 
that 

b{f, D ) = inf b( V ). (B.8) 

If Cd is non-empty, then there exists a unique solution fju given by 

f) D {k) = (M*(k))- 1 ( E X?l?{k)+u D {k)), (B.9) 

0<t<T+T 

where {lf(k))i = E((l t (k))i) and where the multipliers Af, < t < T + T 
and vP(h), < k <f, I < i < N are real numbers such that 

Af>0, (l?,f} D ) H -et>0, *?((l?,i) D ) H -e t )=0, (B.10) 
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for < t < T + T and 

u D >0, f, D >0, uP(k)(f) D (k)) i = 0, (B.ll) 

for < k < T and 1 < i < N. The value of the multipliers X D and v D can 
be determined using the functions F m . 

Next we determine a first approximation and of the multipliers. 

Let 

\V = F-{L{r{u D )-e% (B.12) 

let 

V {1) = B-\ £ tPh + ^l (B.13) 

0<t<T+T 

and let 

» {1) (k)=F Ma{k) ( J2 ^h(k) 

°^ T+T (B.14) 
+ (M a (k) - M b (k))E(f}^(k))) - B{k, l)v {1) {l)), 

for < k < T. The corresponding first approximation fj^ 1 ' of the solution fj 
is then defined by 

°^ f+T (B.15) 
+ (M a (k) - M h (k))E{rp-\k))) -J2B(k,l)f}^(l)), 

l^k 

for < k < T . Formula ( |B.15| ) ensures that r)^ > 0. The approximation 



method defined by formulas (p.!3| ), (P-14 ) and (P-15| ) can in an obvious way 



be generalized to higher order approximations. However the convergence of 
the method must be established. 



C Proofs 

We first give the proofs of the results of §0 and then those of Appendix 
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C.l Proofs of results in §@ 



Proof of Lemma |2.1| Schwarz inequality shows that 77 i— > E(U(oo,r))) is a 
bounded linear map from H to L 2 = L 2 (fi,M), so it is enough to prove the 
statement in the case of a. Since a is maximal, we only have to prove that the 
operator V(a) 3 77 1— > U (00, 77) G L 2 in H, is a closeable, where L 2 = L 2 (Q, M). 
Let {s n } n >i be a sequence in V(a), such that ||s„||-ft — > and \\U(oo,s n ) — 
v \\l 2 0, when n — > 00 for some t> G L 2 . For a G R, let g a be the char- 
acteristic function of the set {u G Q\(J2o<k<f \( u °°(k))( UJ )\ 2 ) 1 ^ 2 ^ a l- Then 
E((g a U(^,s n )) 2 ) < E{a 2 ^<k<fK{k)?) — o 2 ||s n || 2 ^. Hence g a U(oo, s n ) 
converges to in L 2 , when n — ► 00 for every a G R. This gives that v = 0, 
which proves the statement. 



Proof of Lemma |2.2| We prove that (Hi) implies that do is bounded. 
Let X k = rj(k) ■ u°°(k), for < k < f . Then U(oo,rj) = J2 < k < f X k , 

so (E((U(w^)) 2 )) l l 2 < Eo<k<f( E (( x k) 2 )) 1/2 - Since V (k) is ^measur- 
able, since E(\u°°(k)\ 2 \T k ) = E(\u°°(k)\ 2 ) according to hypothesis (Hi) and 
since (X k f) < |M°°(A;)| 2 |r/(A;)| 2 it follows that E((X k ) 2 ) = E(E((X k ) 2 \F k )) < 
E(\u co (k)\ 2 )E(\ri(k)\ 2 ). This shows that a (r]) = E((U(oo,r])) 2 ) < C 2 \\r]\\ 2 H , 
where C 2 = max 0<K f E(\u°°(k)\ 2 ). This proves the statement for o . Since 
< b (7]) = E((U~(oo,r])) 2 ) - (E(U "(00, T]))) 2 < 00(77), it also follows that b 
is bounded. 



Proof of Theorem |2.3| We first prove the statement for o. Because of 
continuity of a, which follows from Schwarz inequality and inequalities ( [2.7|) , 
it is enough to prove inequality ( |2.8| ) for 77 G D(ao) = £f (R , A). In this proof 
we use the notation X k = rj(k) ■ u°°(k), and Y k = J2o<i<k-X-i, for < k < T. 
By the definition of do (see before ( p.6| )) and by the definition of Yf , it follows 
that a (77) = E(Yf). We shall first prove that 

Oofa) > CE( (C.l) 

0<fc<T 

for some C > 0. Let 1 < k < T. For < I < k — 1 we obtain, using that 77(/) 
and rj(k) are T k measurable, using the independence of u°°(l) and u°°(k) and 
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using the T k measurability of E{X k \!F k ) that 

E^Xu) = E(E(( V (l) • u°°(l))( V (k) • u°°(k))\F k )) 
= J2E(m(l)vAk)E(ur(l)u?(k)\F k )) 

= £(£(tj(Z) ■ u°° (1)1^) E(i](k) ■ u°°(k)\F k )) 
= E(E(r)(l) ■ u°°(l)E(r)(k) ■ u°° {k)\T k )\T k )) 
= E(E(rj(l) ■ u°°(l)E(rj(k) ■ u°°{k)\T k )). 

This proves that if < I < k — 1 then 

E{X l X k ) = E{X l E{X k \J" k )). (C.3) 



If 1 < A; < f then E(Y k 2 ) = E((Y k _ x + X k ) 2 ) = E(Y k _ 1 + 2Y k ^X k + X 2 ). 
Formulas (|CJ and E(X 2 ) = E{E{Xl\F k )) now give that 

E(Y 2 ) = E{Yl_ x + TY k -xE{X k \T k ) + E(X 2 k \F k )), (C.4) 

for 1 < k < f . Let a G R N , a ^ 0. Then 

| = |S(o.u°°(A;))| 

i 

< (E(\a ■ u°°(k)\ 2 )) 1/2 = (E(\a ■ u 00 ^ 2 ^) 1 / 2 , 



(C.5) 



where the conditional expectation and the expectation are interchangeable 
because of hypothesis {Hi) and where there is strict inequality since a-u°°(k) 
is a nontrivial random variable according to hypothesis(i?2)- Let 

a(k) = max \E(a ■ u°°(k))\/(E(\a ■ m°°(A;)| 2 ) 1/2 . (C.6) 

|a|=l 

Then < a{k) < 1. Inequality flC5| ) and definition ( |C.6| ) give that 

\E(X k \F k )\<a(k)(E(\X k \ 2 y/ 2 . (C.7) 

It follows from formula (|C.4|) and inequality ( |C.7| ) that 

E(Y 2 ) > E(Y k _ 1 - 2\Y k _ 1 \(E(\X k \ 2 \F k )) 1 / 2 + E{X 2 \F k )) 

= (1 - a(k))E(Y 2 _ l - 2\Y k ^\(E(\X k \ 2 + E(X 2 k \F k )) 
+ a{k)E{{\Y k _ 1 \-{E{\X k \ 2 \F k )) 1 ' 2 ) 2 ) 
>(l-a(k))E(Y 2 _ 1 + E(X^ k )). 
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This proves that there exists C k > 0, (independent of rj) such that 

E(Y 2 )>C k E(Y 2 _ 1 + X 2 k ), (C.8) 

for 1 < k < T. A finite iteration from k = T to k = 1 of this inequality and 
using that Y = X , then proves inequality (|C.1|) . 



Since f]j(k) is T k measurable it follows that 

E ( E **) = E E E ^(^( fc ) E (^(^xwi^)) 

0<fc<T 0<fc<T i,j 

> E J2 E ^(k)v 3 (k)E((u?(k) - E{u?{k))){uf{k) - E(u?(k)))\F k )), 

0<k<T i,j 



where we have used hypothesis (Hi) to deduce the inequality. Then, accord- 
ing to (H 2 ): 

E{ X l) > Y, E CE c ^ k ^)vAk))- (C9) 

0<k<T k>0 i,j 



This inequality and ( |C.1| ) prove that ao(r/) > ^ fc>0 c (k) 2 E(\r](k)\ 2 ) > c 2 ||??|| 2 , 
where c(k) 2 and c 2 are some strictly positive numbers. This proves the in- 



equality for a (rj) in (|2.8j) 



To prove the statement for b, in ( |2.8|) , we first show that if rj e Ti. and 
7] ^ then the random variable U(oo, rf) is not constant (a.e.). Given rj ^ 0, 
let M be the largest integer such that < M < f and r}(M) ^ 0. If M = 0, 
let = JF . If M > 0, let Q be the smallest cx-algebra, such that Tu C 
and such that u°°(k) is ^-measurable for < A; < M — 1. By hypotheses 
and (-££3), and u°°(M) are independent, and X k is ^-measurable for 
< k < M — 1. We remind that, if an absolutely integrable random- variable 
Z is independent of Q, then E(Z\Q) = E(Z). Using hypothesis (.Hi), {H 2 ) 
and (H3) it then follows for the conditional variance that 

EdU^^) - E(U(oo, V )\g)) 2 \g) = E((X M - E(X M \Q)) 2 \g) 

= J2^(M) Vj (M)E((ur(M) - E{u?{M)\g)){uf{M) - E{uT{M)\g))\g) 



hJ 



Yj]i{M)vAM)E{{uT{M) - E{uT{M))){uf{M) - E{u?{M)))) 
J>(M)^(M) % >c 2 \ V (M)\ 2 , 

(CIO) 
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where c 2 > 0. If #(00,77) is a constant, then E((U (00, 77) -E(U(oo, r])\G)) 2 \g) 
= 0, which is in contradiction with inequality (|C.10| ). This proves that 
#(00,77) is not constant, if 77 7^ 0. 

Since #(00,77) is not constant, if 77 7^ 0, it follows that E((U (00, 77) — 
E (# (00, T]))) 2 ) > 0, if 77 ^ 0. This proves that b(rj) > 0, if 77 ^ 0. 

Let c = inf a (^) = i b(rj). Suppose that c = 0. Let {rj n } n >i be a sequence such 
that a(r] n ) = 1 and lim^oo b(V) = 0. Since 6(77) = 0(77) - (E(U (00, t/))) 2 , it 
follows that (E(U (00, i] n ))) 2 — > 1. There is no restriction to suppose that 1 is 
an accumulation point of {i?(#(oo, 77™))}, since the set {r] eH\ a(rj) = 1} is 
invariant under the transformation 77 — > —77. By selecting a subsequence and 
changing the enumeration, we can suppose that E(U(oa,r] n )) — ► 1. Let "H a , 
be the Hilbert space defined on 7i, as a linear space, by the scalar product 
( , ) , where (77, £) = (77, A£)h. The norms a and || ||^ are equivalent, 
due to inequalities (ETTT) and ( |2.8| ) for o and a, so 7i and H a are identical 
as TVS's. Since the unit ball of TC a is weakly compact, (7i is reflexive and 
Banach-Alaoglu theorem, c.f. 0]), it follows that {77 n } n >i has a weakly con- 
vergent subsequence. Once more, by selecting this subsequence and changing 
the enumeration, we can suppose that {77 n } n >i is weakly convergent to an 
element £ G TC. By the definition of weak convergence it then follows that 
1 = lim n ^ E(U (00, 77™)) = E(U (00, £)). This proves that £ 7^ 0. Moreover, 
since £ is an element of the unit ball of TL a and since (E(U(oo, £))) 2 < a (0> 
it follows that a(£) = 1. Therefore b(£) = 0, which is in contradiction with 
the already proved fact that b(£) > 0, for £ 7^ and £ eH. This proves that 
c > 0. Consequently 6(77) > £0(77), for 77 G 7i. The statement for b in (|2.8|) 
now follows from the statement for in 



Proof of Theorem |2.4) Let c = b (resp.a) and C = C\ (resp. C%). C is then 
closed and convex in 7i. Due to Lemma p.2| and Theorem |2]3|, the norm in 
TL is equivalent to the norm given by (c) 1 / 2 . Let TL Z be the Hilbertspace with 
scalar product c. Then C is closed and convex also in TC C - The set C therefore 
has a unique point 77, minimizing its distance (in 7i c ) to the the origin (c.f. 
Theorem 12.3 of §). 

Proof of Lemma |2.5| To prove statement (i), let 77 G Co be a solution of 
equation (0) and let £(#(00,77)) = e. Then E((U (00, fj) -E(U(po, r?))) 2 ) = 
cr 2 , as noted in the paragraph before equation ( |2.4|) . Let inf^ gCl #((#(00,77) 
- E(U (00, 77))) 2 ) = S 2 . Obviously, according to (|2"H), S 2 < a 2 , since 77 G 
C\. Let S 2 < a 2 . Then there exists 77' G C±, such that E((U (00, 77') — 
£(#(00, t/'))) 2 ) = (S') 2 , where (S') 2 < a 2 . Then 77" = (<t/£')V G C and 
E(U (00, 77")) = (er/£')e > e, which is a contradiction since 17 is a solution of 
equation (|1.1|). Therefore cr = S and 77 G Ci is a solution of E((U (00, fj) — 
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E(U(oo,fj))) 2 ) = X 2 . To prove statement (ii), let fj G C\ be a solution of 
equation Q. Let a 2 = E((U(oo, fj) - E(U(po, fj))) 2 ), let m = E(U(po,f))) 
and let M = sup^ gCo E(U(oo, r/)). Obviously m < M according to (|1.1|), 
since ?j G C . Let m < M. Then there exists rf G C , such that cr 2 = 
E((U (oo, i) - E(U(oo, r]'))) 2 ) and E(U(oo, rf)) = M', where m < M' . Since 



e < m < M' it follows that r] 1 G Ci. This shows that equation ( |2.4|) has 
at least two distinct solutions 17 and 77'. This is in contradiction, with the 
fact that the solution of equation ( p^ ) is unique, according to Theorem |2]4. 



Therefore m = M, which proves that fj is a solution of ( |1 . 1|) . 



Proof of Theorem 2.7 Since the two cases are so similar, we only prove the 



statements in case of a. Let C2 be non-empty and let fj be the unique solution 



of equation (|2.12| ). If r\ G C2, t > and 77 ^ 77, then 11(77 + £(77 — 77)) — 0(17) > 0, 
since 17 is unique and C 2 is convex. This gives that (.Da) (17; e) = (A17, e) > 0, 
for e G T C2 (f)), where Tc 2 (fj) is the closed tangent cone of the convex set C 2 at 
17. Let be the set of all e G H, such that (fc) > and (ei(k))(u) = 

a.e. uj G suppfii(k), for < k < T and 1 < 7 < AT. Let if (77) be the closed 
subspace (of 7i) of elements r] e H, such that suppr]i(k) C suppf)i(k), for 
< k < T and 1 < 7 < AT. Then if (77) = ( J + (7y)) , the subset of vectors in 7Y., 
orthogonal to every vector in J+{fj). Let J{fi) be the set of all vectors e 6 W, 
such that (e,r]) n > 0, for all 77 G J+(f])- Then J(r/) = ^+(77) +H(fj). A vector 
e G Tc 2 (fj), if and only if e G J(i7), (m, e)^ = 0, and (lt,e)n > 0, for t G L, 
where L = {t G N | < t < f+T-1, E((AU)(t+l,fj))-c(t)E(K(t,fj)) = 0}. 



Hence 



TcM= f| {e e ^ | (77, e) w > 0}, (C.ll) 

776/(77) 



where I{fj) = J+(fj) U{m, — m} U {i t 1 1 G £}. This is equivalent to 

TcM= D { e ^|(?7,e) w >0}, (C.12) 

776F+(J7) 

where V+(r/) is the closed convex cone generated by I{fj)- This shows that 
Tc 2 (v) = (y+(v))°i the (positive) polar of V + (fj). Since (Af), e) > 0, for all 
e G T C2 (f]), it follows by definition, that A77 G (T C2 (i7))°, the (positive) polar of 
T C2 (fj). But ((^+(77))°)° = V + (fj), because V + (rj) is a closed convex cone. This 
proves that Ar) G V + (fj), so there are /i G R, \t > 0, for t & L, and z/ G J+(fj), 
such that A77 = /jm + XlteL + z/ - Setting A 4 = 0, for < t < T + T — 1 
and t ^ L, we obtain formula ( |2.17| ), since A has an inverse defined on 7i . 



Moreover conditions ( |2.14|) , ( |2.15|) and ( |2.16| ) are also satisfied. 
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C.2 Proofs of results in Appendix 



Proof of Proposition |A.1| The equation (B(k, k) — \)r)(k) = £(k) give, ac- 



cording to formula (|A.2| ), formula (|A.1|) , formula ( |A.3|) and hypothesis (Hi). 
that (B(k,k)-\)r)(k) = (M a (k) - \)r](k) - E(u™ik))(E(u™(k)) ■ E(r](k))) = 
f (Jfe), for < k < N. Since M^(k) = Afg(Jfe) - E{u^(k))E(uf{k)) : it follows 
that (M a (k) - \){r}{k) - E{r](k))) + (M b (k) - \)E{r}{k)) = £(k). Taking the 
expectation of the two members of this expression, we first obtain that 

(M b (k) - X)E(rj(k)) = E{t(k)) 

and then that 

(M a (k) - X)( v (k) - E( v (k))) = £(*) - E(t{k)). 

The matrices M a (k) — A and M b (k) — A are invertible according to the hy- 
pothesis that A i a(M a (k)) and A £ a(M b (k)). The expression flAv§) of 
/c) — A) -1 now follows. 



Proof of Proposition [A.2| Using formula (|A. 12|) it follows that 



9n= 

n+l<r<T 

and using formulas ( |A.9| ) and ( |A.13| ) it follows that 

D n (k) = M a (k)- J2 d r f?N?(k)-\, (C.14) 

n+l<r<T 

where < k < n < T. 

If n = T, then formula (|A.22|) follows from (|A.4|) . Let < n < T and let 



G Hit- The equation B n (k, k)r](k) = £(k) and formula ( |A.20| ) give that 

D n (k) V (k) = ak) + (1 - ■ E( V (k))), (C.15) 

where 7](k) G 7ifc is unknown and < k < n. Taking the expectation on both 
sides we obtain 

D n (k)E(rj(k)) = E(ak)) + (1 - gnMu^ik^Eiu^ik)) ■ E(rj(k))), 



where < k < n. Substitution of ( |C13j ) and ( |C14| ) into this expression gives 
that 

(M a (k) - X)E{r ] (k))-E{u co {k)){E(u QO (k)) ■ E{ V (k))) 

- Yl ^ r 2 (^(fc)E(r / (fc))-E( M 00 (A;))( J E(^(A;))-E(r ? (A;)))) 

n+l<r<T 

(C.16) 
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for < k < n. It follows from ( |A.8| ) that 

(N b n (k)h = (K(k)h - E{uT{k))E{uf{k)), 

which together with QCTTB ) and the definition of M b (k) (see below formula 
(H)) give that 



(M b (k) - d T f r N b r (k)-X)E(r 1 (k)) = E(i(k)), 



(C.17) 



n+l<r<T 



for < k < n. According to hypothesis, A is in the resolventset of M b (k) — 
^ n +i<r<f d rf?N b (k), which proves that 



E( V (k)) = (M b (k) ~ £ drtfN^-Xy^m)), (C18) 

n+l<r<T 

for < k < n. 

Equation ( |C.15|) and the following equation give that 

D n (k)( V (k) - E( V (k))) = - E(t{k)). 

Since A is also in the resolventset of M a {k) — J2 n +i<r<T drfrN°(k), according 
to hypothesis, it follows by expression ( |C.14|) of D n (k) that 



rj(k) - E(rj(k)) = (M a (k) - ^ d r f r N?{k) - A)" 1 ^*) - £(£(*))), 

n+l<r<f 

(C.19) 

for < k < n. Formulas ( |C.18|) and ( |C.19|) prove expression ( |A.22|) . 

Proof of Proposition [A.3| We shall first find a simplified expression of 
certain products, of operators A(k, I) : Hi — > TLk given by (|A.1| ). Let M be a 
real symmetric N x N matrix and let 



I n (k, I) = A(k,n)MA(n,l), 



(C.20) 



for < k < n, < / < n, and < n < T. The operator I n (k, I) : Hi -> H k 
is linear and continuous. Definition (|A.1| ) and formula ( |A.3|) give that 



I n (k,l)r)(l) = E{u™{k){u™{n) ■ (MA(n,l)r)(l)))\F k ) 
= E(u°°(k)(u°°(n) ■ E(Mu°°(n)(u°°(l) ■ V (l))\F n ))\F k ), 



(C.21) 
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for rj(l) e Tii. Since rj(l) is measurable, it follows from hypotheses (Hi) 
and (H 3 ) that 

I n (k,l) V (l) 

= E(u°°(k)(u°°(n) ■ E(Mu^(n)))(E(u 00 (l)\^ n )-r ] (l))\^ k ) 
= J2E(E(u°°(k)(u°°(n) ■ E(Mu 00 (n)))E(ur(l)\^n)Vi(l)\^k) 

i 

= ^(£(u°°(n)) ■ E(Mu°°(n))) 

i 

E(E(u™(k)E(u?(l)\F n )\F n )r,mFk) 
= (E(u°°(n)) ■ S(M«-(n)))^E(E( U -(A;)|^)E( U r(0|^n)^(OI^), 

(C.22) 

for < h < n, < I < n, and < n < f . 

Let k ^ I. The hypothesis (H3) and the first and last member of equality 
([022]) give that 

I n (k,l)ri(l) = (E(u°°(n)) ■ E(Mu°°(n))) 

Y^E(E(u™(k)ur>(l)\F n ) Vi (l)\Fk), k^l. ( C ' 23 ) 

i 

Since r]i(l) is jF n -measurable and n > k, it follows that 

I n (k,l) v (l) = (E(u°°(n)) ■ E(Mu co (n)))A(k,l)r ] (l), k + I, (C.24) 

where < k < n and < I < n. 

Let k = I. The first and last member of equality ( |C.22 ), the ^-fc-measura- 
bility of rj(k) and hypothesis (Hi) give that 

I n (k, k) V (k) = (E(u°°(n)) ■ E(Mu°°(n))) 

j2E(E(u^(k)\^ n )E(ur(k)\^ n )Mk), ( c - 25 ) 

i 

where < k < n. With the notation ( |A.7[ ), formula ( |C.25| ) reads 

I n (k, k)rj(k) = (E(u°°(n)) • E(Mu°°(n)))N^(k)r](k), (C.26) 
where < k < n. 

Next we shall consider products as in (|C.20| ), but for operators B(k,l). 
Let P be a linear continuous operator in 7i„, such that, if r](k) e 7i n and 
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E(rj(n)) = 0, then Prj(n) = Mr](n). We introduce a linear and continuous 
operator J n (k, I) : Hi — » Tik by 

J n {k, I) = B{k, n)PB{n, I), (C.27) 

for 0<k<n, 0<l<n and < n < T. We note that 

E{B(n, l)r](l)) = 0, < I < n. (C.28) 

In fact, since 7/(7) is ^-measurable, it follows from hypotheses (Hi) and (-£^3), 
that 

E(u°°(l) • 77(0) = E(u°°(l)) ■ E(rj(l)), (C.29) 

and then that 

E(B(n, l)rj(l)) = E{E{u™{n){u™{l) ■ V (l) - E(u°°(l)) ■ E^l)))^)) 
= E(E(u°°(n)u? (1)^(1)) - E(u°°(n))(E(u°°(l)) ■ E( V (l))) 

= Y^E{u^{n))E{uT{l))E{r ] m - E(u<* >(n))(E(u<* '(J)) ■ i5fa(0)) 

i 

= 0. 

(C.30) 

If E(rj(n)) = 0, then B(k,n)rj(n) = A(k,n)rj(n), according to (|A.1[) , ( |A.2| ), 
(|C.29|) and hypothesis (-f^i). This gives together with formulas ( |C.27| ) and 
QCg§ that 

J n {k, I) = A{k, n)MB{n, I), (C.31) 

for 0<A;<77, 0</<ri, and < 77 < T. Since definition ( |A.1| ) of B(k, I) 
and (|029D give that 

B(r, s)ri(s) = A(r, 5)77(3) - J E(u 00 (r)) J E(u 00 (s) • r/(s)), (C.32) 

for < r < T and < r < T, it follows from formulas (|C.20| ) and (|C.31|) 
that 

J n (k, l)rj(l) = I n (k, l)n(t) - A{k, n)ME(u x (n))E(u°°(l) ■ 7/(Z)), (C.33) 

for < k < 77, < I < 77 and < n < f. Noting that (use (HJ and (H 3 )) 

A{k,n)ME{u ao {n)) = ^(w 00 ^)^ 00 ^) ■ (ME(u°°(n))))\F k ) 
= E{u°°(k)){E{u°°(ri)) ■ {ME(u°°{n)))), 
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for < k < n, it follows from ( |C.33| ) that 
J n (k,l) V (l) = r{k,l)r]{l) 

-E(u°°(k))(E(u°°(n)) ■ {ME{u°°{n)))){E{u°°{l)) ■ E( V (l))), ( 

for < k < n, < I < n, < n < f and rj(l) G H t . Formulas (|023D and 
(|CT|) give that 

J n (k,l)r]{l) = {E{u°°{n)) ■ {ME{u°°{n)))) 

(A(k, l) V (l) - S(^(A;))(£;(^(0) • E( V (l)))), k + I. 

Definition ( |A.1|) of B(k, I) and and formula ( |C.32|) then give that 

J n (k, l)rj(l) = {E{u°°{n)) ■ {ME{u°°{n))))B{k, l)rj(l), k ^ I, (C.35) 

where < k < n, < I < n, < n < T and r)(l) G Tii. Formulas ( |C.26|) and 
(|CT|) give that 

J n (k, k) V (k) = (E(u°°(n)) ■ (ME(u°°(n)))) 

(N:(k)rj(k)-E(u^(k))(E(u^(k))-E( v (k)))), 1 ' } 

where 0<k<n, 0<n<T and t](k) G Tik- 

Having established ( ggg ) (resp. QCI35D ) and (|C^6|) (resp. ( |C36D ) in the 
case of C = A (resp. B), we now turn to the main part of the proof. 

For < n < T, let P(n) be the statement that the following points (i), 
(ii) and (iii) are true: 

• i) ff, ... , f n , gf, ... ,g n and the matrix elements of D s (k), where < 
k < s and n < s < T, are real valued rational functions of A in IR, 
without singularity in R — cr" +1 . 

If n < T, then Df(T),... ,D n+ i(n + 1) are invertible matrices, for 
A G R — &n +1 , and d , . . . , d n +i are real valued rational functions of A 
in R, without singularity in R — er° +1 . 

(We remind that these functions are defined by (|A.7| ) and (|A.9|) - (|A.13 ) 
and that cr" +1 is defined by ( |A.25| )) 

• ii) If A G R - then C s (k, I), where n<s<f,0<k<f and 
< / < T, given by ( |A.14j )- (|A.20|) , are continuous linear operators 
from Tii to Tik- If moreover n < T, then the operator C s (s, s), where 
n + 1 < s < T, have a continuous inverse in 7i s . 
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• iii) To formulate this part of P(n), we introduce linear continuous oper- 
ators Sf, ■ ■ ■ , S n in H, defined by Sfq = rj and, if n < T and n < s < T, 
by 

(S s r](k)) = v (k), s+l<k<f (C.37) 

and 

(S s 7)(k)) = r)(k) - C s+ \k, s + l)(C s+ \s + l,s+ l))" 1 r/( S + 1), 

(C.38) 

for < k < s + 1. (It follows from statement (ii) of P(n), that the 
operator S s : H — > Ti. is well defined.) 

Let A G R — Then 5 S : 7i — > H. is continuous with continuous 

inverse, for n < s < T, and 

C s = C/ S C (C.39) 



r = ^e, (c.4o) 

for n < s < T, where U s = S s S s+ i ■ ■ ■ Sf. 

We prove, by finite induction, that P(T), . . . ,P(0) are true. The state- 
ment P(T) is trivially true. In fact K — &f +1 — M — = K (see above 
( |A.25| )) and ff, gf, Df(k), where < k < T, given by ( |A.9| ), are well defined 



for A G R. Statement (i) is then true. Since C : — >■ H is continuous, it 
follows that C f (fc, l):Hi^H k ,0<k<f,0<l<f, given by ( fATTj ) and 
( |A.15|) are continuous. Thus, statement (ii) is true. Statement (iii) is also 



true, since Uf = Sf = I, where I is the identity in H. 

Let < n < T. We shall prove that, if P(T), . . . ,P{n) are true, then 
P(n — 1) is true. We remind that a c n+1 C o^j, (see below ([A.26Q ). 

We first prove that (i) of P(n — 1) is true. If n — T, then Df(T), given 
by (|OD , is invertible for A G R-cr£. If n < f, then D T (T), . . . ,D n+1 (n + l) 
are, according to (i) of P(n), invertible matrices, for A G R - tr" C R - 
cr^ +1 and D n {n) is, according to expression ( |C.14j ) of D n {n) and definition 



(|05D of <, also invertible for A G R - <. Thus, D f (T), . . . ,D n (n) are 
invertible matrices, for A G R — er*. The matrix elements of the matrices 
Df(T),... ,D n (n) are, according to (i) of P(n), rational functions (of A), 
without singularities in R — a* (c R — er" + i). Formulas ( |A.10| ) and ( |A.11| ) 



then show that df,... ,d n are rational functions, without singularities in 
R — a".. This proves that the second part of statement (i) of P(n— 1) is true. 
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According to statement (i) of P(n), it follows that ff , ... , f n , gf, ... ,g n 
and the matrix elements of D s (k), where < k < s and n < s < T, 
are real valued rational functions of A in R, without singularity in R — <j" C 
R — er" +1 . Since we already have proved, that d n is a rational function without 
singularity in R — <r", it follows from ( |A.12|) and (|A.13|) , that /„_i, g n _i and 
the matrix elements of P n _i(/c), where < k < n — 1, are rational functions 
of A in R, without singularity in R— <r". This completes the proof of statement 
(i) of P(n- 1). 

Next we prove statement (ii) of P(n — 1). According to (i) of P(n — 
1), which we already have proved, D n (n), is invertible, / n _i, g n -i, are real 
numbers and D n -i(k), where < k < n — 1, are real matrices, for A G 
R - < (c R - Let A G R — <. According to (ii) of P(n), C s (k,l), 

where n<s<T,0<k<T and < Z < T, are continuous linear operators 
from Hi to ft fe . Now it follows from (PO^HPOOQ , that C n ~ l {k, I) : H t -> W A , 
where < k < T and < / < T, are continuous linear operators. This proves 
the first part of statement (ii) of P(n — 1). 

Let A G R — cr^, where c = a (resp. b), when C = A (resp. B). If n = T, 
then a£ = a(M a (f )) (resp. a(M a (T )) U a(M b (f ))) according to (TAT25D 
(resp. ( PT^SD ). We remind that D f (T) = M a (T) - A, (seeJA~|)). According 
to the inversion formula ([OTP (resp. of A f (f ,f) (resp.5 f (f , f )) 

and its hypothesis, it follows that A T (T, T) (resp.-B T (T, T)) has a continuous 
inverse in Tif. If n < T, then (ii) of P(n) gives that C s (s, s) has a continuous 
inverse in Ti. s , for n + 1 < s < T. Since A G R — a^, it follows, in the case of 
C = A from expression (|C.14| ) of D n { n ) an d formula ( |A.21|) and in the case 
of C = B from formula (|A.22|) , that C n (n, n) has a continuous inverse in 7i n . 
This completes the proof of statement (ii) of P(n — 1). 

Finally we prove statement (iii) of P(n — 1). Let A G R — a^. Then A G 
R— 0"n+i) so ^ follows from statement (iii) of -P(n), that 5" s is continuous with 
continuous inverse in 7i, for n < s < T. The operators C n (k,l) : Tii — > Tik, 
< k < T and < / < T and the inverse of C n (n,n) are continuous 
according to (the already proved) statement (ii) of P(n— 1). Therefore, if> n _i 
given by ( |C37[ ) and ( |C38[ ) is a continuous operator in 7i. It is invertible, 
because its inverse is obtained by formulas ( |C.37|) and ( |C.38|) (in the case of 
P(n — 1)) with the sign changed, of the second term of the right hand side 
of ( |C38j ). This proves, that the operators S s : 7i — > H, n — 1 < s < T, are 
continuous with continuous inverse. 

The operator C s : H — > H, satisfy ( |U.39| ), for n < s < T, according 
to statement (iii) of P(n). To establish (|C.39| ) for s = n — 1, it is therefore 
enough to prove that C n ~ l = S n ^C n . Let C = S n ^C n . It follows from 
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QjlgTp and ( P^8| ) that 

C'(k, Z) = C n {k, Z), n < k < f , < I < f , (C.41) 

and 

C'(k,l) = C n (k,l) - C n {k,n){C n (n,n)Y l C n {n,l), (C.42) 

for < k < n and < / < f . To prove that C'(k,l) = C^^M), for 
< k < T and < Z < T, we shall consider (k, Z) chosen according to a 
certain disjoint partition of the set {(fc, I) G N x N | < k < f , < I < f }. 

Formulas ( |OS| ) and ([C^TD show that C'(k, I) = C n -\k,l), for n < k < 
f and < / < f . 

Let < k < n and < I < f . Formula flCggp S iv es that C'(k,n) = 0. 
If Z > n, then statement (hi) of P(n) gives that formula ( |A.17| ), with n + 1 
instead of n, is satisfied, so C n (r, s) = for < r < n + 1, n + 1 < s < T. 
In particular, for the above chosen k and /, C n (k, I) = and C n (n, I) = 0, so 
C(M) = according to ( gH ). Formula ( [ATfD now gives that C"(M) = 
C^-^Jfe, Z), for < < n and < Z < f . 

Since statement (iii) of P(n) is true, it follows from the expression of 
C n (r,s) obtained, by expression ( |A.14| ) in the case when n = T and by 
expression ( |A.18| ), with n replaced by n + 1, in the case when n < T, that 

C n (r,s) = f n C(r,s), (C.43) 

where r^s,0<r<n and < s < n. Formulas (|C.42| ) and (|C.43| ) give that 

C'(k, I) = C n (k, Z) - f n C{k, n)(C n (n, n)y l C(n, Z), (C.44) 

for < k < n and < Z < n. Since A G M — a^, it follows , in the case of 
C = A (resp. B), from formula ( |A.21| ) (resp. ( |A.22| )) and its hypothesis that 
(C n (n, n))^ 1 exists and that (C n (n, n))~ l r](n) = (D n (n))^ 1 r](n), if rj(n) £ H n 
(resp. r](n) G Ti n and E(r](n)) = 0). The product (7(Zv, nXC^n, n)) -1 ^^, Z) 
in (|C.44|) can therefore be simplified, in the case of C = A (resp. B), by 
using the already derived expressions ( |C.24| ) (resp. ( |C.35| )) and ( |C.26| ) (resp. 
([036D ) of J n (/c, Z) (resp. J n (k,l)), defined by ( pjp|) (resp. ( |027|) ), with 
M = (D^n))- 1 (resp. P = (B n (n, n))" 1 and M = (D^n))- 1 ). 

Let k ^ Z, < fc < n and < Z < n. Formulas ( UM) , ( UM) (resp. 
(|C.35| )) and expressions ( |A.10| ) and ( |A.11| ) of d n , give that 

C'(kJ) = C n (k,l)-f^d n C(kJ). 

Formula (|C.43|) gives that 

C'{k,l) = f n {l-f n d n )C{k,l). 
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Expression (|A.12| ) of / n _x and expression (|A.18 ) now show that C'(k,l) = 
C n -\k, I), for k ^ I, < k < n and < I < n. 

Let k = I and < k < n. Since P(n) is true, it follows that C n (k, k) is 
given by ( |A.19|) (resp. ( |A.20| )), with n + 1 instead of n, for C = A (resp. B). 
We first consider the case of C = A. Formulas (|C.44|) , ( |A.19Q and (|C.26|) and 
formulas (|A.10|) and ( A.ll|) of d n , give that 

A'{k,k) = D n {k) - f n d n N«{k)). 

Definitions QQ3|) of D n ^(k) and flP^D of A n_1 (ife, k) now show that A'(Jfe, k) 
= A n ~ 1 (k, k). Next we consider the case of C = B. Formulas flC!13D , flOOl) 
and ( |C36[ ) and formulas (|A.10[ ) and ( |A.11| ) of d n , give that 

B'(k,k) V (k) = (D n (k) - f*d n NZ(k)) V (k) 

- (1 - <7„ - fld n )E{u°°{k)){E{u°°{k)) ■ E( V (k))), 



for n(k) E H k . Definitions (|A~T2l) of g n ^, ([All of D n _i(Jfc) and (|A~20D of 
B n ~ 1 (k, k) now show that B'(k, k) = B n ^ 1 (k, k). This completes the proof of 
formula ( |C.39| ), for n — 1 < s < T. 

The elements £ s G 7i satisfy ( |C39| ), for n < s < T, according to statement 
(iii) of P(n). To establish (|C39 ) for s — n — 1, it is therefore enough to prove 
that C' 1 = S n _i£ n - Let £' = S n _ x £ n . We remind that (C^n,™))" 1 exists 
(see below formula ( |U.44j) ). The definition of 5„_i, by formulas ( |C.37|) and 
(|C.38|) ) (with n — 1 instead of n) gives that 



£'(£)= n<k<f 



and 



e / (fc)=rW-C"(fc,n)(C"(n,n)) 



< fc < n. 



This 



Formulas QOBD and (|A~2^ ) defining f*- 1 , now show that £' = f" -1 . 
completes the proof of statement (iii) of P(n — 1), so P(n — 1) is true. 

The equation C T r\ — £, £ 6 (i.e. equation ( |A.6| )) is, according to 
P(n), (formulas ( |C.39| ) and (|C.40| ) of statement (iii) of P(n)) equivalent to 
C n r] = £ n , < n < T, i.e. equation ( |A.27| ). This completes the proof. 



Proof of Proposition [A.4| Let A e R — o-q and let < n < T. According 
to the definition of cr£ (see ( |AT25|) and ( |AT2oT) ), it follows that A e R - <. 
Expression (|C.14|) of D n (n), formula ( |A.21|) and the definition of <r", show 
that A n (n, n) has a bounded inverse. Formula ( |A.22| ) and the definition of 



cr£, show that B n (n, n) has a bounded inverse. Successive use of (|A.16|) , give 



that C°(k, k) = C k (k, k), for < k < T. Expression ( A. 29 ) now follows from 
( |A.28|) . This completes the proof. 
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